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Abstract 

We extend the lightcone worldsheet lattice approach to string theory, proposed in 
1977 by Giles and me, to allow for coincident D-branes. We find a convenient lattice 
representation of Dirichlet boundary conditions, which the open string coordinates 
transverse to the D-branes satisfy. We then represent the sum over all planar open 
string multi-loop diagrams by introducing an Ising spin system on the worldsheet lattice 
to keep track of the presence or absence of fluctuating boundaries. Finally we discuss 
a simple mean field treatment of the resulting coupled Ising/coordinate worldsheet 
system. The interplay between Neumann and Dirichlet boundary conditions leads to 
a richer phase structure, within this mean field approximation, than that found by 
Orland for the original system with only Neumann conditions. 
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1 Introduction 



The problem of summing planar diagrams is central to many issues in theoretical physics, 
from the large N approximation to QCD [1] to the relationship of string theory to quantum 
field theory as exemphfied by the AdS/CFT correspondence [2]. In the absence of an ana- 
lytic solution to this problem, it would be nice to have an approach amenable to numerical 
simulations. Over 30 years ago, Giles and I proposed a way to digitize interacting bosonic 
string theory, in its first quantized version, as a sum over histories on a lightcone worldsheet 
lattice [3] (GT). Although inspired by perturbation theory, the resulting formalism provides 
a fully non-perturbative dynamics which reproduces the formal^ perturbation theory when 
expanded in powers of the string coupling g. The restriction of this sum over histories to 
planar open string loops should be manageable on a computer. 

In recent years, my colleagues and I have constructed an explicit lightcone worldsheet 
representation of the planar diagrams of a wide range of matrix field theories [6]. These 
constructions rely on fermionic worldsheet ghosts to cancel the bulk degrees of freedom of 
each worldsheet coordinate - the worldsheet systems that reproduce field theory diagrams 
are essentially topological. The ensuing negative signs in the path integrands spell serious 
difficulties for numerical simulations of such worldsheet path integrals. This problem can be 
avoided by replacing each field quantum with a finite tension open string, bringing us back to 
the original GT formalism. Thus we propose that numerical simulations of planar diagram 
sums be performed in the GT formalism at finite string tension, after which conclusions 
about field theory can be drawn by study of the infinite tension limit. 

All open string coordinates in [3] satisfy free end (Neumann) boundary conditions. It has 
long been clear that D-branes [7] provide the key for arranging that the infinite tension limit 
of critical string theory (in 26 or 10 space-time dimensions) yield a quantum field theory in 
lower dimensional space-time. D-brancs arc subspaces on which open strings end, meaning 
that the open string coordinates describing motion perpendicular to the D-branes satisfy 
Dirichlet boundary conditions. The purpose of this article is to explain how such conditions 
can be introduced in the GT worldsheet lattice formalism, and to begin to assess, in the 
context of a simple mean field approximation, their impact on the physics of planar diagram 
summation. 

In addition to bringing in Dirichlet boundaries, the GT lattice formalism also needs 
to be extended to include Grassmann coordinates in order to describe the Ncvcu-Schwarz 
(NS) boson model [8,9], Ramond fermions [10,11], or the superstring [12]. In particular, 
we have suggested that the even G-parity Neveu-Schwarz open string model with SU(A^) 
Chan-Paton factors (NS-I-), which is free of open string tachyons, could be used to establish 
a string representation of large QCD [13,14]. Since the lightcone worldsheet lattice is 
tailor made for summing planar diagrams, its extension to cover the NS-|- open string model 
would provide a promising way to sum planar diagrams on a computer. Then study of the 
To ^ oo hmit should yield new information about large N QCD. 

^Technically Lorentz covariance in the continuum limit requires counterterms to cancel lattice artifacts 
that arise from boundary terms in the integration over moduli (see e.g. [4, 5]). Bulk and boundary terms 
already included in [3] can account for some (perhaps all) of the necessary counterterms. 
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We conclude this introduction with a brief review of the original GT lattice formalism for 
bosonic string theory [3]. It starts with a lattice worldsheet path integral for the lightcone 
quantized free open string [15], 
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where for simplicity we have taken the discrete unit of a to be Tqg (P"*" = MgTq), with a 
the discrete unit oi r {T — {N + l)a). Then a drops out of the formulae and the continuum 
hmit is simply M, N ^ oo with N/M — TqT/P'^ fixed. In the above expression, x1 = xl 
and x^~^^ = xl are fixed by the initial and final states. 

We have written the potential term Vn of the lattice action appropriate to the open string, 
with Neumann boundary conditions which are automatic consequences of the absence of a 
"bond" joining the sites (1, j) to the respective sites (M, j). The closed string action would 
be obtained by simply restoring those bonds. 

v^-'ixi) = vr-'ixD+^f^^^.-xi.r (4) 

More generally, by rearranging the bond patterns in the potential energy term V of the 
worldsheet action S, we can describe any number of closed and open strings. Then interac- 
tions among strings can be achieved by summing over histories in which the bond patterns 
change from time to time [3]. Each appearance or disappearance of a bond is accompanied 
by a factor of g, and each bond interchange by a factor of g"^. 

This general sum over histories would involve wildly nonlocal interactions on the world- 
sheet, which would surely defeat any attempt at numerical evaluation. However, the nu- 
merical prospects are much brighter for the sum over histories corresponding to planar open 
string multiloop diagrams"^. In this case the only changes in bond patterns would be the 
appearance or disappearance of bonds between nearest neighbor sites. Since these changes 
are all local on the worldsheet, techniques of condensed matter physics and quantum field 
theory should apply. For instance, Orland [16] has applied the technique of mean field the- 
ory to study the physics of this planar diagram summation^. He introduced an Ising spin 



^Incorporating an SU(A^) "color" symmetry via Chan-Paton factors, these planar diagrams are singled 
out by 't Hooft's large N limit [1]. With a canonically normalized gauge coupling gg the limit holds ^^iV 
fixed. To simplify writing we will absorb a factor of VN in our coupling: g = gs^fN 

*See [17] for a recent alternative treatment of mean field theory in this context. 
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variable sj — ±1 to represent the two states "on" (s = +1) or "off" (s = —1) of each planar 
bond. Putting P- = (1 + sl)/2, we then have 
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The terms on the last line account for the coupling constant g and the bulk (a) and boundary 
worldshcct countertcrms. Since we have included a fluctuating bond between ?' = 1 and 
i = M, this expression describes the planar evolution of a closed string. For the corresponding 
evolution of an open string, one simply imposes the constraint = 0. 

Although (6) is completely well-defined and finite for fixed M, N, it allows for only the 
two natural countertcrms that were introduced in [3] , and shown to be necessary at the level 
of free strings {g = 0). Indeed, the lattice evaluation predicts that the closed and open free 
bosonic string ground state energies have the M — > oo behavior 
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where G = X^^o(~)"/(2^ + ~ 0.9159656 is Catalan's constant. In [3] we observed that 
the divergent (and non-Lorentz invariant) terms can be absorbed in a and (3. Choosing 

2G{D - 2) 



a — 



+ 0{g% (3 = ^ ln(l + V2) + 0{g') 



(9) 
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gives finite, Lorentz covariant, and correct values for the free string energy spectrum. Al- 
though it was not clearly stated in [3], we must expect that both of these counterterm 
parameters will receive corrections for nonzero coupling g ^ in order to maintain Lorentz 
covariance for D = 26. Indeed, the one loop corrections in bosonic string theory do con- 
tribute to them. An important open problem, not addressed in this article, is to resolve 
whether these two counterterms suffice to render the loop expansion covariant. If not, any 
further counterterms must be identified and incorporated into the formalism. 

The rest of the article is organized as follows. In Section 2 wc present our prescription for 
handling Dirichlet boundary conditions on the lattice. In Section 3, we represent the sum 
over planar diagrams, in which some open string coordinates satisfy Neumann boundary 
conditions and others satisfy Dirichlet boundary conditions, as a sum over Ising spin con- 
figurations, where the Ising spin keeps track of the fiuctuating boundaries. Then in Section 
4, wc give our implementation of the mean field approximation to the Ising spin dynamics. 
Section 5 concludes the article with a discussion of our results and problems for the future. 
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2 Dirichlet Conditions on the Worldsheet Lattice 



For notational clarity we shall use x to describe the string coordinates satisfying Neumann 
conditions, and we shall use y for the string coordinates satisfying Dirichlet conditions. We 
first consider a single free string with Dirichlet conditions — y]^ — 0. The simplest way 
to discretize this string is to use the action: 

N M-l N /M-2 \ 

j=0 1=1 j=l \ i=l / 

Note that 5" involves only M — 1 integration variables for each time shce j. Now consider 
how to obtain this action from the closed string action, which involves M integrations for 
each j, as would be necessary in the sum over histories. Then, the replacement 

{y{ - yir + {y^M - yi-if - (i/if + d/l^-if, (n) 

encounters the problem that the coordinate describes a spurious zero frequency mode. 
An easy fix for this is to include an extra term Tq ^jLid/lf)^ in Sd- 

N M 

^ YEE(^r'-2/i)'+^(^i) 

„ N rM-2 ^ 

Vniyi) - fE E(2/m-2/D' + 2/f + 2/il-i + V^ (12) 
j=i I i=i ) 

Then, since the frequency of this added mode is 0(1) in lattice units, the mode described 
by y'f^ is irrelevant in the continuum limit. But retaining it allows an efficient description 
of the creation and destruction of Dirichlet boundaries without changes in the number of 
degrees of freedom. 

Before turning to that, we give the explicit evaluation of the path history sum for the 
propagation of a free Dirichlet string. Define 

an = 2(A^ + 1) ' n^l,2,...,N (13) 

f3rn = 4sin2— , m = 0,l,...,M-l (14) 

which are the respective eigenvalues of the kinetic and potential bilinear forms occurring 
in 5'n. The eigenvalues of the potential bilinear form appearing in 5*0 are the (3rn, ^ = 
1, . . . , M — 1, plus the eigenvalue 2 for the extra coordinate y-^j^. For economy of writing it 
is convenient to put /3m = 2 and to define the frequencies 

u;^ = 2sinh-^^^, form = l,...,M (15) 
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Then the path integral for a Dirichlet string propagating in d dimensions over time T — 
{N + l)a from y. — to iff — using (12) is 
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where we have taken M large in the last hne and introduced w = e (^+^)'^/^ — e '^^'^o/P'^ 
Here, 
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The first two terms, one linear in M and the other independent of M, contribute divergent 
terms to the continuum P~ = Em /a and violate Lorentz invariance. But they can be 
cancelled by the bulk and boundary counterterms respectively. Here we see explicitly that 
the effect of the mode we added is simply to modify the coefficient (3 of the boundary 
counterterm. Taking the case of D—2 Neumann and 26— D Dirichlet open string coordinates, 
we see that we should have 



2 ~" 1 + a/2 
26 -D 



ln(2 + ^) + 12 ln(l + ^2) + 0(/) 



(19) 



For D — A this \s (3 — —3.91 + 0{g^). Curiously, g — Q, (3 stays negative for < 10 and 
is positive for D > 10. For D ^ W, P ^ 0.041 + 0{g^). 



3 Summing Planar Diagrams 

We turn now to the problem of representing the sum of planar diagrams by introducing the 
same system of Ising spins = ±1 or equivalently P- = (1 + s])/2 = 0, 1 used in the case 
of Neumann conditions [16]. Associate each coordinate yl with the corresponding P- and 
let P- — when the Dirichlet condition apphes to yj. Then we should write the potential 
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term as: 

Y E - y'if + (1 - POPUiy'ii + (1 - PDPi-iuti + 2(1 - if )2/f ] 

-ToY, \yi" - P^PUivl ■ yU ■ (20) 

Then the hghtcone worldsheet action that sums the planar diagrams of Dirichlet open strings 
would be 

N M N M 

^ = YEE(^^^'-2/i)' + ^EE(2/f-2/i-^m^^i) (21) 

j=0 i=l j=l i=l 

Notice that this implementation of Dirichlet conditions has the feature that if every site is 
Dirichlet, i.e. P- — for all then the system is just M independent oscillators with 
frequency of 0(1), and the continuum limit would show no interesting physics. 

In order to describe D = A dimensional physics with a critical string theory in 26 
(bosonic) or 10 (Neveu-Schwarz) space-time dimensions, one can, as in the development 
of the AdS/CFT correspondence [2], introduce a stack of N coincident D3-branes, which 
are 3 + 1 dimensional subspaces on which open strings end. Let us call the 4 coordinates 
parallel to the D3-branes x'^ and the coordinates perpendicular to the D3-branes . For 
the bosonic string / takes 22 values and for the Neveu-Schwarz string it takes on 6 values. 
The coordinates x{a, r) for an open string satisfy Neumann boundary conditions dx/da = 
whereas the coordinates y satisfy Dirichlet boundary conditions y^ — 0. Then a possible 
worldsheet lattice set up would be 

N M 

i,j i,j j=0 i=l 

N M 

+ Y E - ^')'^ + - 22/i ■ vUiP^PU (22) 

j=i i=i 

The sum over planar diagrams is accomplished by summing over all spin configurations 
P- — 0, 1. It is worth pointing out the physical situations represented by the extreme spin 
configurations. If all P- = 1 The x's and y's appear on an equal footing and represent a 
single closed string moving in 25 spatial dimensions (for the bosonic string). In the opposite 
extreme, with all P- — 0, the x's see no potential and represent M free Newtonian particles 
moving in 2 spatial dimensions, and the y's represent M such particles bound by a harmonic 
oscillator potential to the point y = 0. 

The parameters a and f3 characterize the bulk and boundary counterterms respectively. 
They will depend on g( in a way that we do not know a priori. It may well be that coun- 
terterms beyond these will be required to ensure Lorentz invariance, though there remains a 
slender hope that these will suffice. A study of multi-loop corrections in perturbation theory 
will be needed to resolve this issue. As a working hypothesis we shall assume in this article 
that only these counterterms play a role. 
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4 Mean Field Theory 



Mean field theory provides a simple method to understand the physics of large systems, 
although it can be misleading especially near critical points. A convenient framework for 
applying mean field theory to our spin system begins with the addition of a source term 
'Y2ij i^]Pi to the action S. Then, writing the path integral in the presence of k, as Z{k) = 
e--^(«)^ the expectations and correlators of the P's can be obtained as derivatives of J- with 
respect to the k's. Defining 

- iPi) = g. (23) 
it follows that the Legendre transform A{(p) = T — satisfies 

S = -4 

Thus the possible values of are stationary points of the effective action A in the absence 
of sources. 

Up to this point no approximations have been made. The mean field approximation 
consists in replacing the coefficients of the coordinate terms in the action by their expectation 
values: 

pi^m - ^, pipL - {pM-,^ ^ 4>ii (25) 

N M 

S^Y.i^ + - Pi)\ - J2{Pr - if )^ In^ + f E jyxf^' - xif + iy^' - y^] 

i,j i,j j=0 i=l 

N M 

+y E - + - 22/i • yi^M ^ SM (26) 

3=1 i=l 

To do this systematically, one can write S = Sq{<^) + AS{(I)) and treat AS{<p) as a pertur- 
bation. Dropping AS is tantamount to the mean field approximation. From this point of 
view could be chosen to be anything, but it should be chosen to make the perturbative 
corrections as small as possible. Dropping AS* decouples the coordinate path integral from 
the spin sum, so that the whole path integral factors into three parts Z — Z^ZyZgi or 
T — Tx -\- -\- with 

e--^^ = ^M-llif^^^-Pi) + ^pt"-pi?^9+^jPm 

z,{k,) = ^ exp I - ^ (/3(1 - n + (P^+i - In g + k,P^) \ (28) 
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Specializing to static sources Kj — k, the last sum can be thought of as the A'"th power of 
the 2x2 matrix 



T 



e ^ ge 



(29) 



which has eigenvalues 



cosh ^ _ ± \ / sinh^ — + g"^ 



(30) 



The contribution of each eigenstate of T to Zs is weighted by t^, and since — > oo in the 
continuum limit, the + eigenstate will dominate: 



(31) 



Since the coordinate integrations are decoupled from the spin sums at zeroth order in the 
mean field approximation, it follows (for static sources) that 
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Then, making the choices 0j = (P/)o = 4>i and = i^i Pi+i)o = 4^1 optimizes the pertur- 
bation theory in the sense that {AS)q ~ 0. These choices imply a linkage between 0j and Ki 
which can be expressed as 



rvi ^ 



Koi(0) = /3 + 21n 
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(34) 



We have added the subscript to k on the right of the last equation to emphasize that 
that relation between k and (j), which we shall use to eliminate k in favor of in the spin 
part of the effective action, holds only at zeroth order and neglects the back-reaction of the 
coordinate fields-the essence of the mean field approximation. The remaining two terms in 
T are determined at zeroth order by 



J n dxi exp I ^ [(0.^^ - + cp^ixl, - xir] 

ij \ ij 



(35) 
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V, 0.5 




Figure 1: The spin contribution to the effective potential, showing the case of 
g — 1 and P — —1 (lowest curve) /3 = (middle curve) and P — 1 (highest 
curve) . 



Here we restrict attention to the planar (with respect to open string loops) evolution of the 
closed string, so that periodic boundary conditions are appropriate. Then by worldsheet 
translational invariance we can expect that the mean field in the system ground state is 
uniform over the worldsheet (pj — (f). It is then sufficient to take a uniform source: = k. 
We are also interested in the continuum limit M, — > oo, so the dominant contribution to 
the effective action, A = — M{N + 1)k0, will be the bulk term M{N + 1)V proportional 
to the area. The mean field will be determined by minimizing this term, i.e. by minimizing 
the effective potential V(0) = Vg + Vx + Vy. We next list the three contributions to V in the 
mean field approximation, including the —K(f) ~ —K,o{<t^)4> term in the spin contribution^. 



v., = 



•^s(^o(0)) / ,w 1 . / ,x , 

M{N + 1) ^ -lni+-Ko((/>)0 



^ 2^4,(1 - 0) V 40(1 




^The technical details of their derivation can be gleaned from the appendices of [18, 19]. Note that the 
condition Ko{(j)) = determines a stationary point of Vs not V. By minimizing the total V we take the 
back-reaction due to coordinate fluctuations into account in an average way. 
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Figure 2: The spin contribution to the effective potential, showing the case of 
P — and g — 0.1 (highest curve) g = I (middle curve) and g — 10 (lowest 
curve) . 



V, = , ^{D-2) r dxsmh-\^smnx) = V ^"'^"\'f (38) 

M(Ar + l) ^ Vo ^ ^ (2n+l)2 ^ ^ 



(39) 



Note that all of the dependence on the coupling g and the counterterm parameter P is carried 
by Vs. As an ilhistration of this dependence we show a plot of Vs for (? = 1 and three different 
values for f3 in Fig. 1 and for /? = and three different values for g in Fig. 2. Plots of Vx 
and Vy arc shown separately in Fig. 3, and combined (for the case D = 4) in Fig. 4. The 
non-monotonic behavior of this last graph is a direct consequence of the opposite monotonic 
behavior for Dirichlet and Neumann coordinates evident in Fig. 3. These curves will be 
the same for all values of g and (3. Of course the combined plot Fig. 4 will depend on D 
which controls the relative weight of the coordinates with Neumann and Dirichlet boundary 
conditions. Finally in Figs. 5, 6, we plot the total effective potential V = V^; + + for 
the same values of P and g used in Figs. 1, 2. 

In the mean field approximation, the value of controls the effective tension of the string. 
To see how, we write the effective action in the continuum limit: 



Aeff ^ / drda 



x' + y' + TS{cl>x" + ^y') + 4(1 - ^')y' 



(40) 

For oscillations parallel to the D-branes (in the x directions) we have T^J^ — Tq^/^, and 
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Figure 3: The coordinate contributions to the effective potential, Neumann on 
the left and Dirichlet on the right. 



1 



1 




O 0.2 0.4 , 0.6 0.8 1 



Figure 4: The total coordinate contribution to the effective potential, combined 
according to the case D — A. 
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for oscillations perpendicular to the D-branes we have T^^^ = TqcJ). The oscillations perpen- 
dicular to the brane also have an effective mass m^'' ' = v/2(l-02)/a. For generic < 1 
this mass is divergent in the continuum limit and would have the effect of suppressing os- 
cillations perpendicular to the D-branes. For these oscillations to cost finite energy, would 
require (f) — 1 + 0{a^) as a — > 0. 

The plots of the effective potential in Figs. 5, 6 all show that both endpoints = 0, 1 
are local minima. The reason can be seen analytically from the endpoint behavior of V. In 
the — s> limit, Vx approaches zero with infinite positive slope, Vy approaches its limiting 
value with zero slope, and approaches with a negative infinite slope. 

V>^^, V;~0, V:--lni, for 0^0 (41) 

Clearly, dominates as long as D > 2, so V(0) rapidly increases as increases from 0. At 
the other endpoint, = 1, V^^ approaches its value at finite slope, Vy approaches its value at 
negative infinite slope, and Vg approaches its value at positive infinite slope. 

Here Vy, Vs are comparable and we conclude that 

V'~-?^^ln^. (43) 

and V will increase from its value as decreases from 1 as long as D < 26 — 2tt. It follows 
that, for small to moderate coupling, only the endpoints = 0, 1 are candidate minima of 
the effective potential. Which one is actually lower in energy is controlled by the value of 
P, which is not known a priori for general g. From the free open string calculation we know 
that (3 ~ — 3.91 + 0{g'^) for D = 4. For a value of f3 this negative, = is favored over = 1 
by a very wide margin for small to moderate couplings. Keeping /? fixed at this value, one 
would have to go to 51 > 7 for a minimum with < < 1 to develop. Since the "bare" one 
loop correction comes entirely from second order perturbation theory, it should lower the 
zero coupling energy and hence require the O(g^) contribution to f3 to be positive, making it 
a little less negative. In any case, we can safely say that mean field theory predicts that the 
effective string tension will vanish in the system ground state at very weak coupling g <^ 1. 

For g sufficiently large at fixed (3, a local minimum in the effective potential develops 
at some 0o between and 1. However, at the critical coupling where V"(0o) vanishes, 
this local minimum has higher energy than one or both of the two endpoint minima, so 
this minimum initially describes a metastable phase. Near the critical point this metastable 
phase would support finite energy spin waves, signifying the emergence of a new Liouville-like 
degree of freedom on the worldsheet. But eventually for larger coupling, this new minimum 
could become a global minimum and the system ground state would support a finite string 
tension. However, since we do not know P for g — 0{1), we cannot rule out the possibility 
that /? becomes more and more negative as g grows. If this happens the local minimum 
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0.2 0.4 , 0.6 0.8 1 

Figure 5: The effective potential, showing the case oi g = 1 and j3 = —\ 
(lowest curve) (5 — Q (middle curve) and [3=1 (highest curve). 



just described may never become a global one. Thus even in the mean field approximation, 
the jury is still out on the question of whether our lattice model of the bosonic string will 
actually support a finite string tension. 

5 Discussion and Conclusions 

In this article we have shown how to extend the lattice worldsheet formalism for the bosonic 
string to allow for D-branes, and we have applied a mean field approximation to the resulting 
lattice model. Although the mean field analysis allowed us to map out the possibilc phases 
of the system, our a priori ignorance of the value of j3{g) leaves us uncertain about which 
phase is actually realized when g > 0(1). However, the mean field analysis is uneqivocal 
about the weak coupling phase: in it the mean field = and the string tension is quenched 
to zero. 

There are many levels at which our results must be regarded as provisional. First is the 
question: does the lattice formalism accurately represent the bosonic string perturbation 
theory? The answer is yes only if the bulk and boundary counterterms wc have allowed for 
are sufficient to absorb all Lorentz covariance violating artifacts due to the lattice cutoff. The 
evidence for this is so far very meager. These counterterms suffice to render the spectrum 
and tree scattering amplitudes Lorentz covariant. At one loop, we have analyzed the open 
string propagator and shown that for this specific one loop process the counterterms suffice. 
It should not be too difficult to check this conclusion for higher point one loop amplitudes, 
but this has not yet been done. We have no information on this issue at two loops and 
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0.2 0.4 , 0.6 0.8 1 

Figure 6: The effective potential, showing the case of /3 = and gf = 0.1 
(highest curve) g = 1 (middle curve) and g = (lowest curve). 



beyond. Clearly more investigation of this quesiton is called for. 

Even if the worldsheet lattice system is not Lorentz covariant because of the need for more 
complicated counter-terms, it remains a well-defined two dimensional system of scalar fields 
interacting with an Ising spin system. The physics of this system can be analyzed in its own 
right. In this article we have begun this analysis within the mean field approximation. The 
analysis suggests that the system exhibits three distinct phases: a phase with zero effective 
string tension (0 = 0), a phase with maximal effective tension Tq (0 = 1) and a disordered 
spin phase, only stable at sufficiently large coupling, with reduced effective tension Tqv^ 
(0 < < 1). This intermediate phase could potentially support a meaningful infinite tension 
limit in which the effective tension stays finite. In such a limit the planar string diagrams 
would go over into planar quantum field theory diagrams, and one might gain insight into the 
large N limit of certain matrix field theories. All of these conclusions depend on the validity 
of the mean field approximation, which by its very nature is somewhat dubious. But there 
are other approaches to analyzing this system. In particular, Monte Carlo methods seem 
particularly apt, since the path integrand is positive definite and local. Such an analysis to 
test the mean field conclusions would be very welcome. 

Finally, we have to recognize that the tachyon in the bosonic open string theory obscures 
the meaning of the open string loop expansion our model is meant to represent. By itself the 
tachyon could simply mean that the system is being studied in an unstable vacuum, and it 
might disappear once a stable vacuum is found. Since our lattice model is a perfectly well- 
defined physical system, its physics could provide information about the correct stable ground 
state. Indeed, the indication, from the mean field approximation at weak coupling, that the 
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string tension is quenched to zero could be the ultimate fate of the tachyon instability: the 
bosonic open string would then be unstable to decaying into an infinite number of string bits 
[20]. If this is the case, the open bosonic string theory would not be a good starting point 
for understanding large N gauge theory. However, the even G-parity sector of the Neveu- 
Schwarz model is free of open string tachyons and provides a more promising approach to 
large N QCD [13, 14]. This possibility makes the extension of the lightcone lattice worldsheet 
formalism to include Grassmann variables a particularly desirable next step. 
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